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ABSTRACT 

We present observations of the debris disk around 7 Doradus, an FIV star, from the Herschel Key 
Programme DEBRIS (Disc Emission via Bias-free Reconnaissance in the Infrared/ Submilhmetre). 
The disk is well-resolved at 70, 100 and 160 fiui, resolved along its major axis at 250 /zm, detected 
but not resolved at 350 /um, and confused with a background source at 500 /xm. It is one of our 
best resolved targets and we find it to have a radially broad dust distribution. The modeling of the 
resolved images cannot distinguish between two configurations: an arrangement of a warm inner ring 
at several AU (best-fit 4 AU) and a cool outer belt extending from ^55 to 400 AU or an arrangement 
of two cool, narrow rings at ^70 AU and ~190 AU. This suggests that any configuration between 
these two is also possible. Both models have a total fractional luminosity of ~10~^ and are consistent 
with the disk being aligned with the stellar equator. The inner edge of either possible configuration 
suggests that the most likely region to find planets in this system would be within ~55 AU of the 
star. A transient event is not needed to explain the warm dust's fractional luminosity. 
Subject headings: stars: individual (7 Doradus, HD 27290, HIP 19893) - circumstellar matter - 
infrared: stars - submillimeter: stars - techniques: photometric 



1. INTRODUCTION 

Debris disks were first discovered when observations 

with the InfraRcd Astronomical Satellite (IRAS) re- 
vealed that Vega, /3 Pictoris, and Fomalhaut were un- 
expectedly bright at infrared (IR) wavelengths (Aumann 
et al. 1984). Many main-sequence stars have since been 
observed to possess IR emission above the expected pho- 
tospheric level, which is usually attributed to the ther- 
mal emission of dust that is heated by the host star(s). 
The dust is second generation, e.g., produced by ongo- 
ing collisions of larger bodies, since the dust's lifetime 
in orbit is too short for it to be primordial in origin, 
e.g., originating in the protoplanetary disk (Backman & 
Paresce 1993). Models suggest that parent planetesimals 
of at least 10 100 km feed the dust through destructive 
collisions, though how objects are initially "stirred" to 
high enough collision velocities is unclear (e.g., Wyatt 
et al. 2007). One possibility is the formation of bodies 
large enough (> 1000km) to stir planetesimals through 
dynamical interactions (Kenyon & Bromley 2004; Mustill 
& Wyatt 2009). 
Modelling of the stellar spectral energy distribution 
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(SED) at optical wavelengths allows a comparison with 
the observed IR flux, and possible detection of an IR ex- 
cess. Modelling the SED of the excess itself provides a 
measure of the dust temperature and therefore its loca- 
tion (by making assumptions about thc^ emissive proper- 
ties of the dust). This method has proven instrumental 
in relating the properties of debris disks to each other. 
It is generally sufficient to assume that the dust parti- 
cles behave like blackbodies, but this approach does not 
yield the exact disk location because dust at different 
stellocentric distances can have the same temperature 
(grain size and dust location are degenerate) . SED mod- 
eling can also yield incorrect dust properties such as lower 
grain size (compare HD107146, Roccatagliata et al. 2009 
versus Ardila et al. 2004 and Ertel et al. 2011). 

Thus resolving a disk allows for a deeper understand- 
ing of the system since its configuration is directly ob- 
served. Disk sizes determined from resolved images are 
often found to be ^2 5 times larger than those suggested 
by the SED with blackbody assumptions (Schneider et al. 
2006; Wyatt 2008; Matthews et al. 2010; Rodriguez & 
Zuckerman 2012; Booth et al., in press). Resolving the 
disk at more than one wavelength is even more advanta- 
geous, as it can reveal whether the observed configuration 
of the system is wavelength dependent, as it is for Vega 
(Su et al. 2005; Sibthorpe et al. 2010), /3 Leo (Matthews 
et al. 2010; Churcher et al. 2011) and others. For ex- 
ample, a disk may appear larger at longer wavelengths if 
it has two components, and the cooler one (i.e., further 
from the star) dominates at the longer wavelength, as is 
the case for 77 Corvi (Wvatt et al. 2005; Matthews et al. 
2010). 

As more data are collected on debris disks from obser- 
vatories with improved sensitivity and resolution, we are 
moving away from the simple disk models of single nar- 
row dust rings and discovering more complicated shapes 
and configurations. Many debris disks have been found 
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Table 1. Stellar Information for 7 Dor 



P n m f^i" f^T 


Value 




Spectral type 


Fl V 


Gray ct al. (2006) 


R.A. (J2000) 


04:16:01.586 


H0g ct al. (2000) 


Decl. (J2000) 


-51:29:11.933 


H0g et al. (2000) 


PM-R.A. (mas yr"'^) 


101.5 


H0g et al. (2000) 


PM-decl. (mas yi~^) 


184.7 


H0g et al. (2000) 


V magnitude 


4.25 


Balona et al. (1994) 


Distance (pc) 


20.46 ± 0.15 


Phillips et al. (2010)^ 


Age (Gyr) 


0.4'' 


Chen et al. (2006) 


Age (Gyr) 


0.82-2. 19<= 


Vican (2012) 



^Incorporates parallaxes from van Lecuwcn (2007) and van Altcna 
et al. (1995). 

''Estimated uncertainty is a factor of two (0.2-0.8 Gyr) using 
Schaller et al. (1992) isochrones. 

"Estimates of ages determined from chromospheric activity and 
X-ray emission, respectively. 



to contain multiple dust components and/or have ex- 
tended dust distributions with a large range in behaviour 
(e.g., HD 107146: Ertel et al. 2011, /3 Leo: Stock et al. 
2010, C Lep: Moerchen et al. 2007). Some hosts that 
have multi-component disks are also planet hosts (e.g., e 
Eridani: Backman et al. 2009; Reidemeister et al. 2011). 
HR 8799 even plausibly shows the existence of a plane- 
tary system with warm and cold dust that has planets in 
the gap (Su ct al. 2009). This highlights the rich diversity 
of observed planetary systems. 

Since debris disks were discovered with IRAS, sig- 
nificant advancement of the field has been made with 
other observatories such as the Spitzer Space Telescope 
{Spitzer; (Rickc et al. 2005; Bcichman et al. 2006; Bryden 
et al. 2009), the James Clerk Maxwell Telescope (Holland 
et al. 1998; Greaves et al. 1998), and now the Herschel 
Space Observatory (Herschel) by observing the thermal 
emission of the dust. (Optical observations of star light 
scattered by dust are also used to study debris disks.) 
Herschel is well suited for debris disk studies as dust 
emission is well contrasted against stellar emission at its 
wavelength range of 70-500 /im. It is sensitive enough 
to detect the photospheres of nearby stars (and there- 
fore can better determine whether an excess is present), 
and its resolution of 6.7" at 100 fim can probe the sizes 
of nearby disks. We present observations of 7 Doradus 
(7 Dor) and its debris disk that were taken with Her- 
schel as part of the DEBRIS (Disc Emission via Bias-free 
Reconnaissance in the Infrarcd/Submillimetre) Key Pro- 
gramme (Matthews et al. 2010). Stars observed by DE- 
BRIS are taken from the UNS (Unbiased Nearby Stars) 
sample (Phillips et al. 2010), a volume-limited sample, 
and so unbiased toward spectral type, binarity, mctallic- 
ity, and presence of known planets. 

The debris disk around 7 Dor (HD 27290, HIP 19893), 
whose basic parameters are listed in Table 1, was discov- 
ered with IRAS (Aumann 1985; Rhee et al. 2007). It has 
been detected with Spitzer but was not resolved (Chen 
et al. 2006; Koerner et al. 2010). Here we present im- 
ages of the 7 Dor disk taken with Herschel that are the 
first to directly constrain the location of the dust. 7 
Dor is not known to host any extra-solar planets but 
is a target for the exoplanet search using Near-Infrared 
Coronagraphic Imager (NICI) on the Gemini-South 8.1 



m telescope (Chun et al. 2008; Liu ct al. 2010), the results 
of which are not yet publicly available. In Section 2, we 
present the Herschel observations and ancillary data. We 
measure the Herschel fluxes in Section 3. In Section 4, 7 
Dor's photosphere is modeled and used to determine the 
excesses that are observed in the IR and submillimetre 
(submm). We present a basic analysis of the images in 
Section 5 and the more detailed modeling of the SED 
and resolved images in Section 6. The results are dis- 
cussed in Section 7 and we summarize our conclusions in 
Section 8. 

2. HERSCHEL OBSERVATIONS 

Broadband photometric mapping observations were 
performed at 70, 100, 160, 250, 350, and 500 /xm us- 
ing the Photoconductor Array Camera and Spectrome- 
ter (PACS; Poglitsch ct al. 2010) and Spectral and Pho- 
tometric Imaging Receiver (SPIRE; Griffin et al. 2010) 
instruments on board Herschel (Pilbratt et al. 2010). 
These observations were performed as part of the DE- 
BRIS key programme (B. C. Matthews et al., in prepa- 
ration) and were performed in mini and small scan-map 
modes for PACS and SPIRE, respectively. PACS ob- 
serves at 100 and 160 /im simultaneously or at 70 and 
160 /xm simultaneously. Since 7 Dor was observed at 
both 70 and 100 ^m, 7 Dor was observed at 160 /xm twice 
and the map at this wavelength is composed of the data 
from both observations. The coordinates of the peak of 
the emission in the 70 and 100 /xm maps are used to align 
160 /im maps. A summary of the observation parameters 
is given in Table 2. 

The PACS data were reduced using HIPE {Her- 
schel Interactive Processing Environment: Ott 2010) 
version 7.0 Build 1931. The reduction method includes 
some data that were taken while the telescope was re- 
versing the scan direction, and therefore not scanning at 
a constant speed, to decrease the noise in the map. This 
reduces the noise at the map centre (where the target is 
located) by ^25%. PACS data are subject to 1/f noise, 
whose power spectrum is a function of the noise at a 
given angular scale (determined by the scan rate of the 
telescope) . To maximize the signal-to- noise in our maps, 
we filter out signals on angular scales larger than 66, 66, 
and 100" at 70, 100, and 160 /xm in the map generation 
process. The SPIRE data were also reduced using the 
standard Herschel pipeline script in HIPE. 

Figure 1 shows the 60" x 60" region around 7 Dor at 
70, 100, and 160 /im and the 100" x 100" region 
at 250, 350, and 500 /im. A background source, 
which is not listed in the NASA/IPAC Extragalactic 
Database (NED) or the NASA/IPAC Infrared Science 
Archive, is visible ^30" to the southeast. The emis- 
sion from 7 Dor and its disk is well resolved with 
PACS at 70, 100, 160 /im and marginally resolved 
at 250 fxm. The respective full width half maxima 
(FWHMs) of the star-l-disk observations at 70, 100, 
160, 250 and 350 /im are 10.2" x 7.3", 12.7" x 7.9", 
19.2" X 13.1", 26.0" X 19.1" and 27.5" x 22.8" compared 
to the FWHMs for a standard star of 5.6", 6.8", 11.4", 
18.2"and 24.9" (these are the geometric means of the 
beam sizes listed in Table 2). Emission is detected at all 
Herschel wavelengths, however, the background source 
to the southeast, which is well separated from the disk at 
70-160 /um, is harder to separate at SPIRE wavelengths. 
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Fig. 1. — PACS (top row) and SPIRE (bottom row) Herschel observations showing the 60" X 60" and 100" X 100" regions around 7 
Dor, respectively. The beam sizes (the FWHMs of the PSFs) are shown in the bottom right corners. A point-like background source is 
evident ~30" to the southeast of 7 Dor (up is North). The extended emission from 7 Dor is evident in comparison to the background 
source and the beam sizes. The background source is well separated from 7 Dor at PACS wavelengths, but becomes increasingly difficult 
to distinguish from 7 Dor at SPIRE wavelengths. At 500 /im, the flux from both sources is contained within a single beam. 



Table 2. Herschel Observations 



ObsID 


Instrument 


Band 


Time on Target 


Scan Rate 


Date(s) 


Beam Size 


Noise 






(^m) 


i^) 


(" s-l) 


Observed 


{") 


(mjy) 


1342193149/50 


PACS 


100 


1129 


20 


2010 Mar 31 


6.69 X 6.89 


1.3 


1342193149/50 


PACS 


160 


1129 


20 


2010 Mar 31 


10.65 X 12.13 


3.4 


1342220766/67 


PACS 


70 


1129 


20 


2011 Apr 29 


5.46 X 5.76 


1.2 


1342220766/67 


PACS 


160 


1129 


20 


2011 Apr 29 


10.65 X 12.13 


2.7 


1342204956 


SPIRE 


250 


901 


30 


2010 Sep 21 


18.7 X 17.5 


4.9 


1342204956 


SPIRE 


350 


901 


30 


2010 Sep 21 


25.6 X 24.2 


7.3 


1342204956 


SPIRE 


500 


901 


30 


2010 Sep 21 


38.2 X 34.6 


5.3 



Note. — Beam sizes from PACS Observer's Manual v2.3 and SPIRE Observer's Manual v2.4. 



At 250 and 350 /im, the nearby background source and 
7 Dor begin to blend together, whereas at 500 /im, the 
two are indistinguishable. 

3. HERSCHEL FLUX MEASUREMENTS 

Two methods are used to measure fluxes in the Her- 
schel maps: aperture photometry is used when the emis- 
sion is resolved and a point-spread function (PSF) is 
fit to the source when it is unresolved. Thus the flux 
from 7 Dor is measured with aperture photometry at 
70-250 /im, and PSF fitting at 350 and 500 ^lm. The 



southeast background source is consistent with a point 
source in the maps, and so it is fit with a PSF at all 
wavelengths. The PSF fit is done after rotating the in- 
strumental PSF to match the rotation of the telescope at 
the time of the observations^ and using a minimiza- 
tion method. Herschel observations of the diskless star 
a Boo are reduced in the same manner as the DEBRIS 

^ The 160 firci PSF is composed of two equally weighted PSFs, 
each rotated to the corresponding position angle of the telescope 
at the time of the observation. 
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data and used as the instrumental PSF at PACS wave- 
lengths. Empirical SPIRE PSFs are downloaded from 
the ESA ftp site.i*^ 

At 250 and 350 fim, where the background source is 
not well separated from the 7 Dor disk, two PSFs are fit 
simultaneously to the expected positions of each source. 
As a result, the disk flux at 350 fim is very uncertain (as 
well as the background source flux at 250 and 350 /im). 
We quote the flux at 500 /xm from a single PSF fit as 
an upper limit since the background source and the disk 
are within the same beam. The coordinates of the back- 
ground source are determined from the 160 fim map, 
where it is both bright and well separated from the 7 
Dor emission. The position of the PSF fit to 7 Dor is 
fixed to its expected coordinates at the time of the ob- 
servations, given its proper motion (Table 1). This is 
reasonable as there is no evidence of an offset between 
the disk centre and the stellar position given that the 
location of the emission is 1.4 and 2.7" from 7 Dor's ex- 
pected position at 70 and 100 /im (determined from the 
PSF fit) and therefore within Herschel's 2.3" pointing 
uncertainty (Herschel Observers' Manual v4). 

The uncertainty in the flux is determined from the 
noise, listed in Table 2, measured by fitting a PSF to 
400 random locations (the flux is the only free parame- 
ter) within a region devoid of significant emission and of 
similar coverage as the map centre (see B. C. Matthews 
et al., in preparation). (This is scaled by the area of the 
aperture for aperture-measured fluxes.) The PACS cal- 
ibration accuracies of 3%, 3%, and 5% at 70, 100, and 
160 /im, respectively (PACS Observer's Manual v2.3), 
and SPIRE pixel size correction factors and absolute flux 
calibration accuracy of 7% (SPIRE Observer's Manual 
v2.4) are added in quadrature. The measured fluxes of 
7 Dor are listed in Table 3. The fluxes of the south- 
east background source are 11.4 ± 1.7, 20.0 ± 1.9, 35.0 
± 4.6, 28.0 ± 5.7, 18.5 ± 7.9, and <16.7 ± 5.9 mJy at 70, 
100, 160, 250, 350, and 500 /xm. The 350 /xm detection 
is only a 2.5cr detection and the 500 /xm 3a detection 
(with respect to the map noise in Table 2) includes both 
7 Dor and the background source. 

Filtering out signals on large angular scales is neces- 
sary to reduce the noise (see Section 2 for more details), 
however, it also filters out the large angular scales of the 
PACS beam, which extends to about 17 arcmin. As a 
result, the fluxes measured in the PACS maps are too 
low. We correct for this by scaling our fluxes by 1.16 
± 0.05, 1.19 ± 0.05, and 1.21 ± 0.05 at 70, 100, and 
160 /xm (Kennedy et al. 2012). These correction factors 
are determined by comparing the fluxes of bright point 
sources in DEBRIS maps to their predicted photospheric 
flux. It should be noted that these correction factors are 
determined for point sources, but they are reasonable to 
use for an extended source, such as 7 Dor, since the an- 
gular scales of the flltering are still large in comparison 
to the angular scale of 7 Dor's emission. 

4. PHOTOSPHERE AND EXCESSES 

4.1. Modelling the Stellar Photosphere 

Accurate models of the stellar photosphere are crucial 
for debris disk studies, as the analysis is based on emis- 

ftp: / / ftp.sciops.esa.int / pub/hsc- 
calibration/SPIRE/PHOT/Beams / 




0.1 1 10 100 1000 

Wavelength (/xm) 

Fig. 2. — SED for 7 Dor. Herschel fluxes are shown in red. The 
IRS spectrum is plotted in blue. Fluxes which contain a contribu- 
tion from a background source {IRAS fluxes and our 500 fim flux) 
are displayed with open circles. Upper limits are shown with tri- 
angles. The black line traces the photospheric model (Section 4.1). 
The green dotted line traces a sample 70 K blackbody spectrum to 
compare to the shape of the excess spectrum. The solid green line 
traces the total flux from 7 Dor and the sample blackbody. 



sion that is observed in excess of expected photospheric 
emission. The stellar photosphere is modeled using a 
minimization method to flt stellar models from the Gaia 
grid (Brott & Hauschildt 2005) to the observed optical 
and near-IR fluxes (see Table 3). The modeling of the 
stellar photosphere is done consistently for all DEBRIS 
targets (see Kennedy et al. 2012 for more details). 

A photospheric model with a Tcff of 7204 ± 28 K, a 
log(g) of 3.49 ± 0.1, an [M/H] of -0.41 ± 0.17, an of 
1.67 ± 0.22 Rq and an of 6.7 ± 0.1 L© is used. Given 
the uncertainties in deriving log(si) and [M/H] from SED 
modeling, the results are in good agreement (within 0.1 
dex and 0.3 dex, respectively) of Gray et al. (2006). The 
predicted photospheric fluxes from the stellar flt are com- 
pared to the observed fluxes and used to measure the ex- 
cess emission in Table 3. The SED is shown in Figure 2 
with the photospheric model. 

4.2. Spitzer Excesses 

We include Spitzer Infrared Spectrograph (IRS) obser- 
vations of 7 Dor downloaded from the Cornell Atlas of 
Spitzer /IRS Sources (CASSIS). This spectrum is from 
Ardila et al. (2010) and includes AOR 24368640, AOR 
27577600, and AOR 3555584. It is consistent with no 
excess shortward of 22 /im and is therefore in agreement 
with AKARI fluxes that measure no excesses shortward 
of 18 /im. 

5. BASIC IMAGE ANALYSIS 

The size of the disk can be estimated by fltting 2D 
Gaussians to the images. First, a model of the photo- 
sphere is removed from the image by subtracting a PSF 
that is scaled to the expected photospheric emission (Sec- 
tion 4.1). The FWHMs and position angles of the 2D 
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Table 3. Observed Fluxes and Predicted Photospheric Fluxes 
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Wavelength 


Observed Flux 


Instrument 






Method 


Photosphere 


Excess 


Reference 


(fj.m) 


(Jy) 


or Satellite 








(Jy) 


(mJy) 




0.4 


59.8 ± 0.8 


Hipparcos 








59.7 ± 1.1 




H0g et al. (2000) 


0.5 


74.4 ± 0.8 


Hipparcos 








74.8 ± 1.4 




H0g et al. (2000) 


0.6 


71.5 ± 1.4 


Hipparcos 








71.7 ± 1.3 




Ferryman & ESA (1997) 


0.5 


76.5 ± 1.5 


Hipparcos 








75.7 ± 1.4 




Mermilliod (2006) 


1.2 


53.7 ± 12.0 


2MASS 








55.2 ± 1.0 




Cutri et al. (2003) 


9.0 


1.908 ± 0.035 


AKARI 








1.924 ± 0.035 




Ishihara et al. (2010) 


12.0 


1.202 ± 0.051 


IRAS 








1.179 ± 0.022 




Moshir et al. (1990) 




(mjy) 










(mJy) 






18.0 


514.2 ± 21.0 


AKARI 








494.1 ± 9.1 




Ishihara et al. (2010) 


23.7 


315.6 ± 3.2 


MIPS 






PSF fit 


286.0 ± 5.3 


29.6 ± 6.2 


K. Su, private communication 


25.0 


292.3 ± 22.0 


IRAS 








256.7 ± 4.7 




Moshir et al. (1990) 


60.0 


196.7 ± 11.0 


IRAS 








44.2 ± 0.8 




Moshir et al. (1990) 


71.4 


170.7 ± 8.1 


MIPS 






PSF fit 


31.1 ± 0.6 


139.6 ± 8.1 


K. Su, private communication 


70.0 


171.0 ± 8.7 


PACS 


30 


X 


19" aperture 


31.1 ± 0.6 


139.5 ± 8.8 


This work 


100.0 


<476.7 


IRAS 








15.7 ± 0.3 




Moshir et al. (1990) 


100.0 


148.4 ± 7.7 


PACS 


30 


X 


17" aperture 


15.7 ± 0.3 


132.7 ± 7.8 


This work 


160.0 


134.3 ± 14.1 


PACS 


30 


X 


20" aperture 


6.4 ± 0.1 


127.9 ± 14.1 


This work 


250.0 


52.5 ± 6.5 


SPIRE 


30 


X 


24" aperture 


2.45 ± 0.05 


50.0 ± 6.5 


This work 


350.0 


23.5 ± 8.0 


SPIRE 






PSF fit 


1.24 ± 0.02 


22.2 ± 8.0 


This work 


500.0 


<16.7 ± 5.9^ 


SPIRE 






PSF fit 


0.60 ± 0.01 


<16.1 ± 5.9* 


This work 



^Although there is a 3a detection at 7 Dor's position at 500 /xm, the flux is listed as un upper limit since there is a known background 
source within the 500 //m beam. 



Table 4. 2D Gaussian Disk Fits 



Band 
(//m) 


Beam'' 
(") 


^obs.maj 


^obs.min 

('^) 


^act,maj 
(") 


^act,min 
{") 


Inclination 

(°) 


Position Angle 

n 


70 


5.61 


12.7 ± 0.7 


7.9 ± 0.4 


11.4 ± 0.8 


5.6 ± 0.6 


61 ± 6 


56 ± 6 


100 


6.79 


14.2 ± 0.7 


8.2 ± 0.4 


12.5 ± 0.8 


4.6 ± 0.7 


68 ± 5 


52 ± 4 


160 


11.36 


20.2 ± 0.8 


13.3 ± 0.5 


16.6 ± 1.0 


6.9 ± 1.0 


65 ± 5 


62 ± 5 


250 


18.2 


26.3 ± 4.2 


19.3 ± 3.1 


19.0 ± 5.8 


b 


>52'' 


55 ± 22 


350 


24.9 


26.3 ± 6.9 


22.7 ± 6.8 


c 


c 


c 


c 



^Tlie effective beam size (geometric mean) is listed and used to deconvolve the sizes using Equation (1). 

''The disk is only marginally resolved along one axis at 250 /im. Therefore Sact.min cannot be derived for the short axis and consequently 
only the upper limit on the inclination can be calculated. 
■^The system is consistent with being unresolved and symmetric at 350 /im. 



Gaussian models are listed in Table 4. The quoted error 
on a fitted parameter gives the range of parameter val- 
ues for which the value is within 10% of the minimum 
value. The value is measured in a small region 
around 7 Dor.-'^^ 

The major and minor FWHMs are deconvolved from 
the beam size using 

^obs = ^act + ^beam) (1) 

where (?act is the actual angular size of the object, ^^bcam 
is the effective beam size cLiid. ^obs is the observed angu- 
lar size of the object. (This is the relation for convolving 
a Gaussian with a Gaussian.) Assuming the true shape 
of the disk is azimuthally symmetric, the deconvolved 
major and minor FWHMs are used to estimate an incli- 
nation of ~65° from a face-on orientation (see Table 4) 
that is consistent with the more detailed image modeling 
(see Section 6.3). 

The fitted parameters of the 2D Gaussian models 
are used to define elliptical annuli to measure the sur- 

" 30"x30", 30"x30", 40" x 40", 70" x 70", and 
80" X 80"regions are used at 70, 100, 160, 250, and 350 /im, 
respectively. 



face brightness profiles of the disk in the 70, 100, and 
160 fj,m images where the disk is resolved along both the 
major and minor axes. Figure 3 shows that the shape of 
the profiles is wavelength dependent. This is consistent 
with a broad radial distribution of material as dust fur- 
ther from the star will be cooler and dominate the flux 
at longer wavelengths. This is also reflected in the larger 
FWHMs at longer wavelengths listed in Table 4. The 
underlying dust distribution is investigated by modeling 
the images in the following section. 

6. DISK MODELING 

6.1. Basic Model 

We first present the underlying set up of the debris disk 
modeling that we implement in Sections 6.2 and 6.3, to 
introduce the various parameters. The simple approach 
to modeling debris disks is to assume that the dust is 
contained within a narrow/discrete ring at some distance 
from the host star. This model is easily extended to 
model multiple narrow/discrete rings by summing up the 
individual contributions of each ring using a modified 
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Semimajor axis (arcsec) 

Fig. 3. — Surface brightness profiles of tfie 7 Dor dislt. The 
average surface brightness in eUiptical annuli with increasing semi- 
major axis at 70 (blue triangles), 100 (green squares) and 160 (red 
circles) fira. la errors are shown. The stellar model is subtracted 
from the maps to measure the profile of the disk (solid lines). (The 
model of the nearby background source is also subtracted from the 
maps.) Beam sizes (dotted lines) are plotted for comparison. The 
shape of the profile is consistent with a radial distribution of dust 
where the cooler dust is brighter and the warmer dust is fainter at 
longer wavelengths (see the text). The surface brightness profiles 
fall off to similar levels beyond ~10" at all wavelengths. 

blackbody function (Dent et al. 2000), 

= 2.35 X 10"" d-^ ct{R) B^{\,T{R)) X^^ (2) 

R 

where F^^ is the flux density (in Jy) at wavelength A 
(in /im), d is the distance to the star (in pc), and <t{R) is 
the cross-sectional area (in AU^) of the ring at radius R 
(in AU). B^{\,T{R)) is the Planck function (in Jy sr"!) 
for the ring at radius R with temperature T (in K). 
Equation (2) is modified by X\, where X\ — (A/Aq)'^ 
for wavelengths longer than Aq and X\ — 1 otherwise. 
This accounts for the fall off of the spectrum that is 
observed to be steeper than the blackbody function at 
submm wavelengths for most debris disks (Dent et al. 
2000). 

If the dust grains act like blackbodies, the radius of a 
dust ring can be derived from its temperature given the 
luminosity of the star, L^: 

Rbb = 278.3^ lV^ (3) 

However, using the above method to model the SEDs of 
debris disks typically underestimates the observed sizes 
of the disks by a factor of 2-5 (Schneider et al. 2006; 
Wyatt 2008; Matthews et al. 2010; Rodriguez & Zucker- 
man 2012; Booth et al., accepted, and references therein) 
because grains with a lower radiation efficiency, depen- 
dent on their size and composition, will be hotter than 
blackbody grains at the same distance would. 

Equation (2) can be expanded to model dust in a 
wide/extended belt by treating the belt as a series of 
dust rings and setting a and T to follow power-law dis- 
tributions. The cross-sectional area is determined from 



the optical depth, t{R), given by 

t{R) = To {R/R^)-\ (4) 

where tq is the optical depth at Rq = 1 AU and 7 de- 
scribes the fall off of the surface density at further dis- 
tances from the star. Similarly, T(i?) is given by, 

r(i?) = To (R/Ro)-', (5) 

where Tq is the temperature at i?o = 1 AU and S de- 
scribes how quickly the temperature of the dust grains 
declines at larger radii from the star. The blackbody 
temperature distribution given in Equation (3) has the 
values of (5 = 0.5 and Tq = 438.4 K for 7 Dor. The 
cross-sectional area, a, the fractional luminosity, /d, 
and the optical depth, r, of the dust are related by 
a{R) = 2Tr RdRT{R) = /d47ri?2 (see Krivov 2010 for 
a review). For a narrow/discrete ring model, we assume 
dR/R^ 0.1. 

This construct is valid for modeling both the SED and 
the resolved images. The SED modeling traces the tem- 
peratures of the dust within the disk (e.g., T, a), whereas 
the image modeling constrains its physical arrangement 
(i? and therefore the dependence of T and a on R). This 
approach (used for all DEBRIS modeling papers) param- 
eterizes the shape of the SED at each radius. How that 
shape relates to the physical parameters of dust compo- 
sition and size distribution will be discussed in a later 
paper. 

6.2. SED Modelling 

For the majority of debris disks, only unresolved pho- 
tometry is available. We include an analysis of the SED 
without the resolved spatial information to highlight how 
the resolved images reduce the ambiguities. A single dust 
temperature is unable to account for the spectral breadth 
of the excess emission (Figure 2). We therefore model the 
SED with two different extended temperature distribu- 
tions: a model of two discrete dust rings and a model of 
a single extended dust belt. All models have the same 
fractional luminosity of 2.6 x 10~^. 

The SED model of two discrete rings accounts for the 
broad shape of the SED with a warm ring (T ^ 125 K, 
Rbb ^ 12 AU) that dominates the flux at IRS wavelengths 
and a cool ring (T ~ 50 K, Rbb ^ 77 AU) that accounts 
for the Herschel fluxes. This model implies that the Her- 
schel images would be dominated by dust at a single 
temperature. 

The extended belt SED model requires a relatively 
flat proflle of the cross-sectional area (7 ^ in Equa- 
tion (4)) to account for the similar fluxes at 70, 100, and 
160 ^m and the mid-IR (IRS and MIPS 24) excesses. 
Using a blackbody temperature distribution, this model 
suggests that the dust extends from ~13 AU to ~166 AU. 

The immediate contrast between unresolved and re- 
solved photometry is evident in the predictions from the 
SED models. The SED model of two discrete rings sug- 
gests that the Herschel fluxes are dominated by a single 
narrow ring at 77 AU, however, the resolved images show 
an extended distribution of dust that cannot be modeled 
with a single ring, as shown in Figure 4. The discrep- 
ancies between what would be interpreted from the SED 
and what is revealed by modeling resolved images (see 
the following section) are clear. Resolved imaging is the 
only way to constrain the possible spatial distributions. 
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Fig. 4. — An attempted single narrow ring model of the 7 Dor disk. Left to right: the observations, the simulated Herschel observations 
of the model, the residuals, and simulated observations of the model at high resolution at 70 iim. The color scale of the residuals extends 
from — 3fT to 3cr. Clearly the single narrow ring model cannot adequately explain the observed emission and we must turn to a model that 
includes dust at a range of radii. 



Table 5. 



Fitted Parameters for the Two Narrow Rings 
Model of the Herschel Images 



Table 6. Fitted Parameters for the Wide Belt Model 
OF the Herschel Images 



Parameter 


Value 


Ri 


70 AU 


R2 


190 AU 


T2-1 


0.85 


Inclination 


71° 


Position angle 


55.9° N of E 


70 /im flux 


140 mjy 


100 iJ.m flux 


142 mJy 


160 /im flux 


125 mJy 


6.3. Image Modelling 



To keep the results from SED modeling and image 
modeling distinct, we describe our imaging models as 
"narrow rings" and "wide belt" which are analogous to 
the "discrete rings" and "extended belt" SED models, 
but of course with different fit parameters. We also re- 
fer to the dust observable to Herschel as "cool dust" 
and dust that is responsible for any excesses at shorter 
wavelengths (IRS and MIPS 24) that does not effect the 
Herschel observations as "warm dust." 

The extended spatial distribution of dust is clear in 
the images. The best-fit imaging model of a single nar- 
row dust ring is unable to account for the emission on 
the largest and smallest scales (see Figure 4). Although 
the images show that all the dust cannot lie at the same 
stellocentric radius, the resolution of Herschel is not suf- 
ficient to distinguish between different extended config- 
urations. Below, we show that both a model of two nar- 
row rings of dust and a model of a wide belt of dust 
are able to reproduce the observations of the cool dust 
component. Both models are parameterized by an inner 
radius and an outer radius. (A constant opening angle 
of 10° was used and reasonable variations to this do not 
affect the fitting.) The narrow rings model describes two 
discrete dust components at these radii with no dust be- 
tween them whereas the wide belt model is described by 
a smooth distribution of material between these edges. 
We do not intend to fully investigate the possible pa- 



Parameter 


Value 


Inner radius 


55 AU 


Outer radius 


400 AU 


To 


3.4 X 10-5 


7 


0.8 


To 


585 K 




0.5 


Inclination 


68.5° 


Position angle 


55.6° N of E 



^These parameters are fixed in the model. 



rameter space given the limited resolution. Rather, we 
compare the two models to demonstrate the uncertain- 
ties on the radial extent of the disk. These two extended 
spatial distributions are reasonable to consider as they 
are observed in other resolved debris disks (e.g., 77 Corvi: 
Wyatt et al. 2005; Matthews et al. 2010, HR 8799: Su 
et al. 2009, HD 181327: Schneider ct al. 2006). The two 
configurations are modeled with similar techniques at 70, 
100, and 160 /im, where the 7 Dor disk is well separated 
from the nearby background source. A gi'id analysis is 
used for the two narrow rings model (as in Booth et al., 
accepted), and a combination of by-hand and min- 
imization (as in Kennedy et al. 2012 and Wyatt et al. 
2012) is used for the wide belt model. 

The best-fit imaging model of two narrow rings is de- 
termined using a grid of parameters: the radius of the 
inner ring (i?i), the radius of the outer ring (-R2), the in- 
clination of the disk, and the ratio of the optical depths 
of the two rings (t2_i = ti/t2, where ti and T2 are the 
r for the inner and outer rings, respectively).^^ The grid 
size (and therefore computation time) is minimized by 
fixing the position angle of the disk to that determined 
from the 2D Gaussian model of the 70 fim image (Sec- 
tion 4). Additionally, for each set of grid parameters, 
the total disk flux, and therefore the total r of the sys- 

12 The grid contains 225,000 models with radii from 55 to 330 AU 
and T2_i from 0.1 to 1.5. 
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Fig. 5. — Imaging models of the 7 Dor disk, [top to bottom) The observations, the simulated Herschel observations of the model, the 
residuals, and the simulated observations of the model at high resolution at {left to right) 70, 100 and 160 /im for the imaging models 
of the cool dust. The narrow rings model is on the left, and the wide belt model is on the right (see the text). The color scale for the 
residuals map extends from — 3(t to +3cr. The residuals show that both models of the cool dust successfully reproduce the resolved PACS 
observations. 



tern, is fit using the package MPFIT (Markwardt 2009). 
The model and its residuals in Figure 5 show that the 
two narrow rings model is able to reproduce the Her- 
schel observations. The fitted parameters for this model 
are listed in Table 5. 

The wide belt imaging model is parameterized by an 
inner radius, a disk width, a power-law surface density 
with slope described by 7 (Equation (4)), an inclination 
to the line of sight and a position angle on the sky. The 
temperature of the dust grains is modeled to fall off with 
6 = 0.5 (Equation (5)) as it was not necessary to deviate 
from this power law. The images are modeled using a 
Levenberg-Marquardt minimization. This technique 
runs the risk of finding a local minimum in the value, 
however, it is less computationally intensive than using 
a grid and the best-fit model reproduces the images, and 
so must be considered a plausible representation of the 
disk structure. The model and its residuals are shown in 
Figure 5 and the fitted parameters are listed in Table 6. 
The fitted Tq (585 K) is hotter than that for a blackbody 



(438 K) . This is likely due to the dust grain composition 
and/or size distribution. 

We compute the reduced of each model in the 
49" X 49", 49" X 49", and 50" x 50" region" around 7 
Dor at 70, 100, and 160 iim. The residuals of both 
models are reasonably low with Xdof ~ ^.78 and 0.71 
for the narrow rings and wide belt models, respectively. 
This suggests that both configurations are possible for 
the cool dust around 7 Dor. We use the best fits of each 
model to estimate the uncertainty on the radial extent of 
the disk. The narrow rings model shows the significant 
emission extends out to at least 190 AU. Emission be- 
yond that is predicted by the wide belt model to extend 
out to 400 AU, however it is difficult to constrain given 
its low surface brightness. This makes the outer radius 
of the wide belt model highly uncertain. There is also a 
discrepancy of approximately 15 AU between the inner 

The region around the background source is not included in 
this calculation as each model treats the background source differ- 
ently. 
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Fig. 6. — SED of the 7 Dor disk derived from the imaging mod- 
els. Photometric excesses are displayed with black circles and IRS 
excesses are in grey. The narrow rings imaging model is plotted in 
red and the wide belt imaging model is plotted in blue. Individual 
components of each model are plotted with dashed lines and the 
solid lines trace the total flux of each disk model. Both models re- 
quire two components to reproduce the observations. The narrow 
rings model agrees well with the MIPS 24 fim data and places the 
rings at ~70 and ~190 AU. The wide belt imaging model that ex- 
tends from ~55 AU to '~400 AU cannot account for the excesses at 
shorter wavelengths. Therefore, an additional component, a warm, 
narrow ring with a temperature of ^225 it 100 K, is needed and 
its properties are estimated from the unresolved mid-IR fluxes. 



extent of each model. This demonstrates that the dust 
visible to Herschel has an inner radius around 60 AU, 
but cannot be constrained to better than 25%. The res- 
olution of this disk is sufficient to place much tighter 
constraints on the geometrical viewing parameters of the 
disk. The inclinations of the imaging models (68.5° and 
71°) are in good agreement and both models measure a 
position angle of 56°. 

Figure 6 shows the SEDs generated for the imaging 
models. We have adopted the standard /3 = 1 that is 
found for most debris disks (Dent et al. 2000) and fixed 
Ao = 160 /im to accommodate the SPIRE observations. 
Both the narrow rings and the wide belt imaging mod- 
els overpredict the observed 100 /^m flux. This discrep- 
ancy is also present in the SED models (Section 6.2) and 
therefore not an issue with the imaging models them- 
selves. The apparent flatness of the SED between 70 and 
160 /im is not characteristic of blackbody emission. More 
sophisticated models could explore whether this feature 
is due to the dust's composition, causing broad emission 
features from minerals or ices, or due to the shape of the 
size distribution. For example, Lebreton et al. (2012) 
show different SED shapes for HD 181327's debris disk 
for various grain compositions and size distributions. 

The SED of the wide belt imaging model suggests that 
the excess at shorter IRS and MIPS 24 wavelengths is due 
to a separate component of warm dust that does not con- 
tribute to the emission at Herschel wavelengths, rather 
than an extension of the belt. We therefore attribute the 
remaining warm excess in the context of this model to 
a separate warm component of dust and estimate it to 
have a temperature of 225 K and a fractional luminosity 



of 4.4 X 10^^. In order to minimize the warm compo- 
nent's contribution to the flux at Herschel wavelengths, 
we adopt a Ao of 40 ^m. The parameters of this com- 
ponent are highly uncertain given the low significance of 
excess and the fluctuation of the degree of excess between 
different IRS reductions, and so we estimate an uncer- 
tainty of ^100 K on the dust temperature. Given that 
the excess is apparent in all available reductions of the 
IRS data and MIPS photometry at 24 /im, which is even 
more discrepant, we are confldent that it exists. The 
225 K dust temperature suggests that the warm inner 
ring lies at ~4 AU under blackbody assumptions using 
Equation (3). The SED modeling supports that such a 
warm dust temperature is present even though its physi- 
cal location is poorly constrained and could lie anywhere 
between 2 and 12 AU given the uncertainty in tempera- 
ture. Typically, the blackbody radius underestimates the 
true size of the disk by a factor of up to ~5 (Schneider 
et al. 2006; Wyatt 2008; Matthews et al. 2010; Rodriguez 
& Zuckerman 2012; Booth et al., accepted) implying the 
true radius could exceed the nominal blackbody radius 
even further. 

Modelling the SED and resolved images suggests that 
there are at least two dust components in 7 Dor's debris 
disk. However, the resolution of the images and the SEDs 
are unable to distinguish between a configuration of two 
narrow dust rings (at 70 and 190 AU) and a configuration 
of a wide outer dust belt (from 54 to 400 AU) and a warm 
inner ring (around 4 AU). 

7. DISCUSSION 

7.1. Properties of ^ Dor's Disk 

It is important to constrain where the dust is to de- 
termine the implications for a possible planetary system 
(Section 7.2). The modeling of the resolved images shows 
that the dust population probed at Herschel wavelengths 
could be distributed in either a wide belt and a narrow 
ring or two narrow rings. Figure 6 shows that both mod- 
els have two dust components. The narrow rings model 
requires a gap between the star and the inner narrow 
ring at ^70 AU and another gap between the inner ring 
and the outer narrow ring (at ^^190 AU). Whereas the 
wide belt model requires a gap between a narrow ring at 
~4 and the wide belt that begins at ~55 AU. All models 
have an estimated fractional luminosity of ~10~^, within 
the upper limits that are typically found for other debris 
disks around F-type stars (Moor et al. 2011). 

A common clearing mechanism that is proposed for 
such gaps, for example, is the presence of planets, how- 
ever, they would be hard to form at such large distances 
from the star (e.g., Weidenschilling 2003). Therefore, the 
wide belt is physically more plausible since it does not 
require planets at such large radii. However, there are 
problems with the wide belt model as well, as it would 
require a very broad disk of dust producing planetesi- 
mals, rather than containing them in narrow rings. 

The inclination of the disk is ^69° from a face-on ori- 
entation. This suggests that it may be aligned with the 
stellar equator given the inclination of 7 Dor itself is 
^70° (Balona et al. 1996), though the orientation of the 
stellar inclination is not known. 

The fractional luminosity of the warm inner ring of 
~4.4 X 10^^ is about the maximum that can be expected 
to arise from coUisional processes for Sun-like stars from 
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Wyatt et al. (2007) 

/max = 0.16 X 10-3 R'/^ M"'/'' L-^'^ t-^^, (6) 

where R is the dust radius in AU, tage is the age of the 

star in Myr, is the mass of the star in solar masses, 
and is the luminosity of the star in solar luminosities 
(Section 4.1). Given 7 Dor's age of 400 Myr, /max is 
2.5 X 10"^ for a dust ring at 4 AU and therefore transient 
events are not needed to explain this amount of hot dust. 

7.2. Relation to Planetary Systems 

7 Dor is an example of how we are continuing to un- 
cover more structure in debris disks with the resolution 
and sensitivity of modern observatories to move away 
from the simple pictures of single narrow dust rings. Ob- 
serving 7 Dor with Herschel demonstrates the impor- 
tance of surveys across different wavelength regimes. In 
the case of 7 Dor, the Spitzer and Herschel observations 
considered independently would give an incomplete pic- 
ture of the system. Our understanding of the debris disk 
around 7 Dor now is that it is more complex than a sin- 
gle ring. We find that it is composed of dust at multiple 
radii by showing two extreme scenarios of this are both 
possible: two cool, narrow rings beyond 70 AU or a cool, 
wide belt beyond 55 AU and a separate warm inner ring. 
In both cases the dust is located between ~55 AU and a 
few hundred AU. Depending on the exact scenario (par- 
ticularly in the case of a wide belt), an additional warm 
component of dust (at several AU from the star) may 
be required to reproduce the mid-IR excess. The wide 
belt model with a warm inner component is much like 
the Solar System, although it is much younger than the 
Solar System, the dust is much brighter and no planets 
in the system have been detected yet. 

7 Dor's multicomponent disk has implications for the 
possibility of the system hosting planets. For planetary 
systems, such as the Solar System and HR 8799, that 
host both planets and a debris disk, the planets are ob- 
served to lie within or between the regions of dust (Moro- 
Martm et al. 2010). This suggests that gaps in debris 
disks are good places to look for planets, particularly if 
the dust ring has a sharp inner edge (e.g., Fomalhaut; 
Kalas et al. 2005, 2008). Resonance overlap studies have 
shown how a planet can sculpt the inner edge of a disk 
(e.g., Quillen 2006, Mustill & Wyatt 2012). Migrating 
planets as well as planets with highly eccentric orbits 
can create gaps and asymmetries in debris disks (see 
e.g., Wyatt 2003). 7 Dor is a good target for the At- 
acama Large Millimeter/Submillimeter Array (ALMA) 
as 7 Dor's distance is appropriate to resolve a sharp in- 
ner edge. Such observations would constrain the location 
of edges and therefore help distinguish between the two 
models presented here. If 7 Dor's disk is indeed best de- 
scribed by a wide belt, the best place to look for planets 
would be within ~55 AU of the star, where the gap be- 
tween the inner ring and the wide belt is predicted. Sim- 
ilarly in the context of the narrow rings model a planet 
could be sculpting the inner ring at 70 AU and/or planets 
could be responsible for the gap between 70 and 190 AU. 

Determining the inner edge of the disk could also 
have implications for the histories of the planetary sys- 
tems. For example, Wyatt (2003) present a model of a 
Neptune-mass planet migrating from 40 to 65 AU. Dur- 
ing the migration, the inner edge of the disk is pushed 



outward and objects in the disk are swept into resonances 
with the planet. Migration of the Solar System plan- 
ets has been proposed as a possible trigger for the Late 
Heavy Bombardment (LHB) , an epoch of intense crater- 
ing responsible for many of the features we still see on the 
Moon today. The LHB severely depleted the Solar sys- 
tem's debris disk rendering such a system unobservable 
at the distance of 7 Dor within the current detection 
limits and available observatories (Booth et al. 2009). 
Although it would be unlikely to observe a planetary 
system while its planets are migrating, because of the 
short timescales over which this takes place (on the or- 
der of tens of Myr) , the inner edges of debris disks and 
the amount of dust within them are tools to consider the 
Solar System's history in the context of other planetary 
systems. For example, the fact that there is so much dust 
still around 7 Dor suggests that it has not undergone an 
instability as destructive as the LHB. 

8. SUMMARY AND CONCLUSIONS 

We observed the debris disk around 7 Dor with Her- 
schel and detect emission at all six wavelengths of PACS 
and SPIRE. The disk is well resolved at 70, 100, and 
160 /um. It is resolved along its major axis at 250 fiui; 
at 350 fxm, the emission is point-like. The emission at 
500 lira cannot be separated from that of a nearby back- 
ground source. Our measurements arc consistent with 
the disk being aligned with the stellar equator given the 
inclination of ~70° that is measured for 7 Dor based on 
its stellar oscillation modes (Balona et al. 1996). 

The SED of the dust emission has a shape that is too 
broad to arise from dust at a single temperature. There is 
both cool dust (observable to Herschel) and warm dust 
(evident from IRS and MIPS 24 excesses). Similarly, 
the resolved images cannot be modeled by a single nar- 
row dust ring, suggesting any temperature distribution 
within the cool dust arises from a radially broad distribu- 
tion of dust and not dust at a single stellocentric radius 
with a distribution of grain properties. 

We have modeled the resolved images, and therefore 
cool dust, at 70, 100, and 160 /xm as arising from two 
narrow dust rings and as a single wide dust belt. The 
narrow rings model has two rings at large distances (70 
and 190 AU) from the star and accounts for all the ob- 
served excess across mid-IR and submm wavelengths. 
The wide belt model has a cool, wide belt (extending 
from 55 to 400 AU) and accounts for far-IR and submm 
excesses, but not IRS and MIPS 24 excesses. Conse- 
quently, a warm inner ring must be included in this model 
to account for the flxixes detected at shorter wavelengths. 
Therefore, both models require a total of two dust com- 
ponents (with a total /cj ^ 10"^) in the debris disk: two 
narrow rings of cool dust, or a wide belt of cool dust and 
a narrow ring of warm dust. As both models produce 
reasonably low residuals, the Herschel resolution is un- 
able to distinguish which configuration best represents 
the 7 Dor disk. 

In the context of the wide belt imaging model, the IRS 
and MIPS 24 excesses are attributed to a warm inner 
ring. Given the variability between the conclusions of 
different IRS observations and reductions, we do not at- 
tempt to place strong constraints on the properties of this 
ring, but estimate its temperature to be ^^225 ± 100 K 
and its fractional luminosity to be -^4.4 x 10"^. We de- 
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rive a blackbody radius of ^4 AU for the warm dust ring 
given a 225 K temperature, although it will be larger in 
the case of non-blackbody grains. The fractional lumi- 
nosity of this ring is within the levels that are expected 
for the steady state evolution of a ring at 4 AU and there- 
fore a transient event is not needed to explain the levels 
of warm dust in this system. 

Planets arc observed to lie between dust components in 
systems where both dust and planets have been observed. 
7 Dor is therefore a good candidate for planet searches, 
particularly if planets are responsible for the lack of dust 
between the inner and outer dust components. The most 
likely region to find planets would be within 55 AU of the 
star, which both models support as clear of dust. Planets 
found beyond 55 AU would offer support to the narrow 



rings model as they could be responsible for the lack of 
dust between the two rings of cool dust. 
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